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ABSTRACT 

For main-sequence stars beyond spectral type M5 the characteristics of magnetic 
activity common to warmer solar-like stars change into the brown-dwarf domain: the 
surface magnetic field becomes more dipolar and the evolution of the field patterns slows, 
the photospheric plasma is increasingly neutral and decoupled from the magnetic field, 
chromospheric and coronal emissions weaken markedly, and the efficiency of rotational 
braking rapidly decreases. Yet, radio emission persists, and has been argued to be 
dominated by electron-cyclotron maser emission instead of the gyrosynchrotron emission 
from warmer stars. These properties may signal a transition in the stellar extended 
atmosphere. Stars warmer than about M5 have a solar-like corona and wind-sustained 
heliosphere in which the atmospheric activity is powered by convective motions that 
move the magnetic field. Stars cooler than early-L, in contrast, may have a jovian- 
like rotation-dominated magnetosphere powered by the star's rotation in a scaled-up 
analog of the magnetospheres of Jupiter and Saturn. A dimensional scaling relationship 
for rotation-dominated magnetospheres by Fan et al. (1982) is consistent with this 
hypothesis. 

Subject headings: stars: magnetic fields — stars: low-mass, brown dwarfs — stars: 
late-type — planets and satellites: general 



Introduction 



Main sequence stars have a convective envelope around a radiative core from about spectral 
type A2 (with an effective temperature of T e jf 9000 K) to about M3 (T e jf 3200 K). Stars 
warmer than about F2 (T e fj ~ 7000 K) have but a shallow convective envelope that sustains at 
most weak magnetic activity. As one of these stars, the Sun exhibits many of the properties of 
the magnetically-driven variability of the population of F2-M3 dwarf stars. The radiative losses 
from their chromospheres, transition regions and coronae generally decrease with age as stars lose 
angular momentum through a magnetized wind. The radiative losses from these distinct thermal 
domains scale through power laws with the average magnetic flux density on the stellar surface. 
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All of these stars exhibit signatures of flaring that increase with increasing quiescent activity. The 
radio emission from all of these stars, albeit generally weak, scales lineary with the coronal X- 
ray emission, sugg estive of persistent non-ther mal energetic-particle popul ations within their o uter 
atmospheres (e.g.. ISchrijver and Zwaanl . l2000l . and references therein; see iBerger et all l2008al . for 
a recent version of the radio/X-ray scaling). 



Stars cooler than about M2 to M4 are expected to be fully convective (jChabrier and Baraffd . 

2000), lacking the convective overshoot layer and tachocline that are thought to be important 
to the solar dynamo. Yet these stars are capable of generating magnetic field, with a rotation 



actiy i ty relationship tha t cont inues clearly to at least spectral type M8 (e.g., iMohantv et al. 
20021 . iReiners and Basril . 120091 ). However, their magnetic field appears to be a predominantly 
axi-symmetrical lar ge-sc ale poloidal field with a t most a slow evolution in the surface pattern 
(jDonati et all 120081 : also lReiners and Basril . 120091 ) . in contrast to warmer stars that show relatively- 
rapidly evolving, non-axisymmetric fields evolving subject to flux emergence and differential rota- 
tion. 

Spectra of FeH lines reveal that stars at least down to M9 have magnetic fields, and that 
the coolest stars beyond M6 for which fields can be measured hay e magnet ic flux densities of 
fB = 1.5 kG or more (e.g.. IBerger et all l2008bl . lHallinan et all 120061 and 2007. IReiners and B asri. 



1999 



2007. X-ray flaring has been reported down to at least M9.5 (e.g., 



Liebert et al. 



199 



Fleming et all 1200.1 IBerger et all |2008a|) . 



9. 



Reid et al. 



Beyond about M8, the relationship between chromospheric (Ha) or coronal ( X-ray) emission 
and rotation rate b eco mes a weak tendency and ey en that disappears beyond L0 (e.g. . IReiners and Basri 
2009, iBasril . |2009j; see IMohantv and Basril . 120031 . for a discussion of difficulties in assigning spec- 
tral types). Whereas essentially all stars of spectral type M8 have strong H a emission, this drops 



to 60% at LP, k, 15% at L4, and to less than than 10% by L5 (e.g., Mohantv et al. 



Mohantv and Basri 
2007h . 



2003 



West et al 



2004 



Reiners 



2007, 



Reiners and Basril 



2002 



Schmidt et al. 



In contrast to the weakening of the rotation-activity relationship for the traditional chromo- 
spheric (Ha) and coronal (X-ray) indicators of magnetic activity beyond about M8, the radio 



lumin osity continues to increase with increasing angular velocity up to at least LP (jBerger et al. 



2008al ). with (often variable) radio emission detected to at least L3.5 (|Bergerl . I2PP6I ). Interestingly, 
the essentially linear relationship between radio and X-ray luminosities for all other cool stars 
breaks down beyond M5, with the X-ray luminosity dropping by a factor o f about 3, POO below that 
relationship at a given radio luminosity by spectral type L0 (IBerger et all l2008al ). 



The reduction of the traditional diagnostics for chromospheric and coronal activity, the weak- 
ening of their dependence on rotation rate despite the presence of strong magnetic fields, and the 
increase in time scale for rotational braking (discussed in Sect. 4), may have their origin in the 
very low degree of ionization of the photospheric plasma. IMohantv et al\ (120021 ) argue that beyond 
M5, the photospheric plasma is weakly coupled to the magnetic field, and the high resistivity and 
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associated diffusivity of the plasma should render it very difficult to generate electrical currents by 
con vection-driven field motions or to transmit any generated in the stellar interior to the corona (see 
also iMohantv and Basril . 120031 . who also discuss the role of dust in these cool atmospheres) . The 
persistence of radio emission into the L-type range, despite these changes in activity characteristic 
of all other cool stars, suggests we explore a fundamental change in the character of stellar magnetic 
activity: the weakening of a well-developed convection-powered stellar corona and associated stellar 
wind may cause the stellar outer atmosphere to change from a solar-like corona with heliosphere 
to a Jupiter-like magnetosphere shaped by the incoming wind of the interstellar medium (ISM) 
impinging on the stellar magnetic field, while powered by the stellar rotation. 



2. Radiative and rotational energy losses 

Coronal emission and magnetic braking both weaken rapidly beyond mid-M. Up to about spec- 
tral type LO, the energy lost by rotational braking appears to be consistent with that expected from 
a solar-like coronal domain beneath an outflowing wind, as can be seen from the following argu- 
ment. I start from the premise that the stellar atmosphere is essentially hydrostatically stratified. 
Near the stellar surface, an isothermal plasma has a pressure scale height 

akTm p Rl 

Hp = GM* ' (1) 

with stellar radius R*, temperature T, proton mass m p , and gravitational constant G. Here, I 
assume a pure hydrogen gas for an order-of-magnitude estimate. The constant a = 1 for a neutral 
gas, and a = 2 for fully-ionized hydrogen. A largely-neutral photospheric gas of T = 2000 K on a 
compact ultra-cool star with surface gravity g ~ 10 5 ' 4 ~ 9gg cm/s 2 would have H p (2000 K) ~ 9 km. 

In order to lead to magnetic braking, plasma of sufficient density needs to exist out to the 
lesser of the Alfven radius (to which corotation of star and high-atmospheric plasma is enforced) 
and the radius of synchronous rotation (beyond which centrifugal forces dominate). The latter 
(which is likely the smaller of the two for rapidly-rotating L-type dwarf stars) is 

*-(^)», 

for stellar mass M* and angular velocity O = 2tt/P, at rotation period P. The characteristic value 
of Rs ~ (see Table 1) is so much larger than H p (2000 K) that the plasma density at Rs would 
be insufficient to lead to significant angular momentum loss. 

For a plasma at coronal temperatures the density scale height is much higher, of course. Only 
very few stars beyond M9 have been detected in X-rays (e.g.. lAudard et all 120071 . iRobrade and Schmittl , 



2008). Their quiescent emissions are at a level of about Lx/L^ Q \ ~ 10 4 , or log(Lx) ~ 25.5, while 



others have even lower upper limits. 
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Assuming a hydrostatically stratified atmosphere, the electron density n e associated with such 
an emission can be estimated from the X-ray luminosity Lx based on the plasma's volumetric 
emission, n^A(Tc), and characteristic volume, AirRlH p (Tc), 

L x = ^RlH p (T c )n 2 e A{T c ) (3) 

for a plasma at temperature Tc, with an emissivity A(Tc) ~ 2 x 10~ 18 T -2 / 3 (for log (Tc) G 
[5.5,7.5]); for Tc = 1.5 MK, comparable to that of the bulk of the solar wi nd and consistent with 
the v ery-limited X-ray spectral infor mation available on L - type dwarf stars (jRobrade and Schmitti . 



20081 ). A(1.5MK) ~ 10~ 218 (e.g., ISchriiver and Zwaanl . |2000| ). The coronal-base density thus 
estimated for L-type stars with log(Lx) ~ 25.5 and Tc = 1.5 MK is uq ~ 10 8 ' 8 cm~ 3 . 

To estimate the plasma density at height Rs within the equatorial plane, the density for an 
isothermal atmosphere can be approximated from a balance between pressure gradient, gravity, 
and centrifugal acceleration: 



dn(r) _ / GM*m p \ n(r) fQ 2 m p \ 



dr 



V 2kT c ) + \2kT c ) 



+ [— r - !i )n(r)r (4) 



(note that the s tratification of a subsonic wind is very similar to a static stratified atmosphere, e.g., 
Hansteenl . 120091 ) . With the above value of no, n(Rs) ~ iq 5 - 6-5 - 8 cm -3 for the characteristic L2 and 



L5 stars in Table 1 (the centrifugal force below Rs modifies the density profile by a factor of < 2 
relative to the stratification in the absence of rotation). 

Now assume that the centrifugal force beyond Rs accelerates plasma outward, which is replen- 
ished from below at (at most) the thermal velocity 

( kTc \ l (*\ 
vth = 7 , (5) 



for adiabatic index 7, so that the upper limit for the mass loss is 

M = a47rRlv th m p n(R s ), (6) 

with a a geometry factor. If the mass loss were isotropic, a = 1, but as mass is probably lost from 
only part of the surface area, I use a = 1/2 below. 

The stellar magnetic field can enforce corotation out to a distance Ra where the plasma (3 
becomes of order unity, given by 



B { R*\ 



U\R { J MRA)kT c , (7) 

if the field is approximated by a dipole of characteristic strength Bq at the stellar surface; n(r) is 
assumed to be given by Eq. (j3J). 
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Outflowing plasma in a spherical shell carries an angular momentum 




(8) 




(9) 



With the values in Table 1, one finds Wc ~ 10 25 ' 8 erg/s for log(Lx) ~ 25.5. For the characteristic 
L2 star in Table 1, this compares well to the estimated loss W ro t from rotation- age studies (see 
Sect. 4), so that we can conclude that the existence of a faint thermal corona in early L-type stars 
is compatible with the inferred magnetic braking. 

The estimatd Wc is while about an order of magnitude larger than W ro t for the M7 dwarf 
star. This may be because of the exponential dependence Wc on Tc, or because the corona of the 
M7 star may be solar-like, in the sense that a few relatively bright and dense regions dominate the 
emission, so that the base density from Eq. ([3]) would be an overestimate for the highest coronal 
loops. 

The estimates in this section obviously depend on the assumed values of Lx and Tc- The value 
of Wc scales with Lx about as a square root, so that our conclusions are not strongly modified 
if the characteristic value of Lx observed for the L-dwarf binary Kelu-1 is somewhat larger than 
what is characteristic of the ensemble of similar L-type dwarf stars. The dependence on Tc is much 
stronger, not surprisingly, because of the near-exponential stratification of density in units of the 
pressure scale height: for T c = 1 - 3MK, log(W rot ) = 24.2 - 27.3. 



Even though magnetic braking remains compatible with the concept of a hot corona, the radio 
emission in ultracool dwarfs is disproportionally strong. One possible cause for this that should be 
explored, and can then be ruled out, is the interaction with the ISM. Observations suggest that 
there is relatively little heated coronal plasma well into the L-type spectral range, there should 
not be much of an associated stellar wind either. The motion of a star with at most a weak wind 
relative to the ISM should lead to the formation of a magnetopause. The energy input into the 
stellar magnetosphere from this interaction is expected to be low, based on the following. 

The standoff (or Chapman-Ferraro) distance of the magnetopause is set by where the dynamic 
pressure p w v ^ of the ISM wind with density p w and relative velocity v w equals the magnetic pressure 
of the stellar field. Observations suggest that the large-scale field of very cool dwarf stars can be 
approximated by a dipole field, so that the field scales with p/r 3 , with the magnetic moment 
{j, = B^Rl, where is the field strength near the stellar equator. Rcf for a late M-type or L-type 
dwarf star can thus be estimated from p w v^ oc B 2 (Rcf)/8tt, or: 



3. The magnetosphere of the coolest dwarf stars 



Rcf 



3.5( 



B% 



g q/6 



(10) 



R* 



n v 100 
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The constant of proportiona lity is determined by using the observation that for Earth Rcf,® ~ 
10 Earth radii (jRusselj . 120071 ). for a wind sp eed of v mo ~ 4 (in uni ts of 100 km/s) and a characteristic 



solar-wind particle density np = 10 cm 3 dFeldman et al 
polar field strength ([Allenl . 11972 ). 



19771 ). and B G ~ 0.6 G for the Earth's 



The relative motion of stars through the ISM average s about 40 km/s for a sample of very cool 
dwarf stars (e.g., IWood et all 120051 . ISchmidt et all 120071 ) . For a hypothetical dwarf star that has 
no stellar w i nd, th at moves through an ISM like that around the solar system with no = 0.1 (e.g., 
Wood et all |2002| ) at i>ioo = 0.4, and Bq = 1.5 kG, one finds Rqf,*/R* ~ 32; that relative scale 
is similar to the geometry in the case of Jupiter's magnetosphere, for which Rmp/Rj = 42 (e.g., 
Walker and Russelll . 119950 • 



To estimate how much power the ISM wind might impart onto a stellar bowshock and poten- 
tially into the stellar magnetosphere, we estimate the power in the bulk kinetic energy over the 
cross section of the bow shock with cross section ttR cf : 



1 n 20 3 ( R CF \ ' 

10 u n vf o 



R* 
R Q 



(11) 



so that at uq = 0.1, wioo = 0.4, and R* = O.3i? , P w m 10 20 ~ 22 erg/s for characteristic M7 to L5 
stars. This lies orders of magnitude below the Ha and X-ray luminosities (see Table 1) and below 
the power that needs to be extracted from the stellar rotation with age (see below). Hence, the 
ISM wind is not a significant source of magnetospheric activity in dwarf stars through mid-L, even 
if it shapes a close-in asteropause. 



4. From a stellar corona to a planetary magnetosphere 

As the solar-like corona fades away from mid-M to mid-L, one may expect increasing signa- 
tures of a rotation-dominated magnetosphere towards late L-type stars like that for Jupiter. In 
such a magnetosphere, energy is taken from star's rotational energy through a torque applied by 
outflowing plasma . The energy loss, W TO t, can be estimated from the rotation-age relationship. 



Reiners and Basril (|2008l ) infer time scales for magnetic braking that increase from about 3 Gyr for 
an M7 star to 7 Gyr for an L2 star, and over 10 Gyr for an L5 star (the latter is relatively poorly 
constrained and, perhaps, much larger). Using 

W, M = Ai/fi* (12) 

yeidls W ro t ~ 10 26 ergs/s (Table 1). In the early L-type range, this is comparable to Wc from 
Eq. ([9]) that could be taken away by the outflow from the relatively weak corona, i.e., that would 
be plausibly consistent with the existence of such a weak corona. But both W ro t and Wc exceed 
X-ray and Ha losses, i.e., the rotational energy loss exceeds the radiative losses, in contrast to what 
is seen in warmer solar- type stars. 
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Let us now explore an qrder-of-m a gnitu de scaling for the power expected from a rotationally- 
dominated magnetosphere. iFan et al\ (|1982l ) argue for a scaling based on a dimensional analysis 
('principle of similitude') that relates the magnetic moment, Mb = B£ 3 , the rotation period, P, 
and the radius, R a , at which the acceleration process begins to the power, W , generated: 



w 



Att 2 K M| 
c 2 P 2 ii a ' 



(13) 



where the constant K has to follow from a measurement until a full theory is developed. Based 
on the properties in Table [H using a mean photosperic field strength of 1.5 kG, and with Rg for 
R a , values are found fo r 10 log (if ■ Wmt/ W) of 0.6, 0.4, and 0.6 for the characteristic M7, L2, and 
L5 stars, respectively. IFan et al\ ( 1982 ) estimate a range of values for the constant 10 log(K) for 
Jupiter from —1.7 to —0.7. The value of K for Jupiter lies, remarkably, within the range of values 
needed to let W ro t and W be comparable for the stars in Table 1. In contrast, this scaling applied 
to the Sun yields log(V^Q) ~ 21.3 — 22.3, which is well over a thousand times less. At least W for 
ultracool dwarfs is much larger than the solar value, while with the only available calibration point 
- Jupiter - one concludes that W is close to both W ro t and Wc for an early L-type star. 



5. Discussion and conclusions 

Somewhere along the spectrum of stars, brown dwarfs, and planets a transition from an outflow- 
driven asterosphere to a field-shielded rotating magnetosphere must occur. Based on the observa- 
tional evidence, I argue that this transition occurs at the bottom end of the true stellar range of the 
main sequence: for stars cooler than about spectral type M5, the properties characteristic of solar- 
like activity progressively disappear, while, in contrast to warmer stars, magnetic braking extracts 
more energy from the star than needed to power the chromospheric Ha and (weak) coronal X-ray 
emissions. Objects cooler than about L0-L2 may exhibit outer-atmospheric phenomena similar to 
those in the rotation-driven magnetosphere of Jupiter rather than to those in convection-driven 
corona and wind of the Sun. 

The energetics of the main auroral oval of Jupiter dominate over the phenomena associated 
with the lower-latitude structures that are connected directly with the movement of Io, Europa, 
and Ganymede through the jovian magnetic field, and those associated with the higher-latitude 
polar-cap emissions that appear to be driven by the solar wind. Jupiter's main auroral oval is not 
- in contrast to Earth's - the separator between open and closed planetary magnetic field in the 
interaction with the interplanetary magnetic field, bu t instead is mapped to lower-latitude closed 



magnetic field (see, e.g., the summarizing discussion by lCowley et q/1 12003 ). This auroral structure 
is thought to be generated by the precipitation of energetic electrons created by the electric current 
system that is involved in a phenomenon referred to as "corotation" (or rather the breakdown 
thereof) in the middle magnetosphere: beyond the interface where the jovian magnetic field is 
strong enough to enforce corotation of the plasma, a current-system is induced that includes a 



- 8 - 



disk-shaped near-equatorial extrusion in which field and plasma interact through Lorentz forces to 
extract energy from the planet's rotation. 

If the rotation of the coolest dwarf stars similarly powers their magnetosphere, this energy is 
ultimately drawn from the energy of rotation. In the case of the Jupiter, Io and Ganymede provide 
plasma conveniently high in the magnetosphere, outside the distance of geosynchronous rotation 
but within the Alfven radius. In the case of the ultra-cool dwarf stars, it is of course possible 
that one or more close-in planets act as a similar plasma source, but rather than postulating such 
planets, I hypothesize that plasma is provided by a tenuous hot stellar corona formed either by 
residual solar-like activity (perhaps associated with the overturning field of a turbulent dynamo) or 
magnetospheric activity (somehow formed by the breakdown of corotation). How, quantitatively, 
this balance changes between LO and L5 remains to be established. 

For stars up to late-M, the power extracted from the stellar rotation lies below the characteristic 
Ha and X-ray emissions. The estimated rate W ro t at which rotational energy is lost from an L2-type 
dwarf, in contrast, is larger than the characteristic outer-atmospheric losses b ased on the coronal 



X-ray luminosity. In the scenario of the rotationally-dominated magnetosphere, Eviatar and Siscoe 



(|1980l ) argued that at most half of the total power taken from Jupiter's rotational energy would 
be available for radiative emissions, including the aurorae. The rough estimate of the radiative 
losses of L-type dwarfs (as reflected by Lx) is not inconsistent with that argument, but further 
observations are needed to test this in more detail. 



The quiescent radio emission from F-M type stars appears to be predom i nantl y 



chrotr on emission from high-energy coronal electrons (see, e.g., the review bv lGiidell . 12002 ) . lHallinan et al. 



gyrosyn- 



(|2008l ) argued that the radio signal of several very cool dwarf stars (M8.5, M9, and L3.5) is associ- 
ated with the electron cyclotron maser (ECM) process; they base that argument on the observation 
that the radio signal is periodically fully circularly polarized, whereas there is also an unpolarized 
component which they attribute to depolarized ECM emission because its brightness temperature 
is incompatible with incoherent synchrotron radiation. They pointed out that the ECM mecha- 



nism is thought to be responsible for kHz-MHz emission from Jupiter (cf., the review by IZarka 



19981 ). Thus, both the enhanced radio to X-ray ratio and the radio polarization signature support 
the analogy of the magnetospheric processes around Jupiter and around ultracool radidly-rotating 
dwarf stars. 

The plasma transport in the magnetosphere of Jupiter, is thought to occur by convective 
cells triggered by the centrifugal inte rchange instability, in which cool, den se volumes change places 



with hotter, less dense volumes (e.g. JHill et al.l . l2005l . lChen and Hill 120081 ) . In the case of ultracool 



dwarf stars, the stellar equivalent of this process into the stellar magnetotail shaped by the ISM 
wind could be the process by which angular momentum is removed from the star. 

In view of the arguments discussed in this paper, it appears warranted to consider the possibility 
of a very similar system of currents in the case of main-sequence stars cooler than about LO as in 
the jovian magnetosphere. Perhaps much, if not most, of the Ha emission and the relatively weak 
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coronal X-ray in stars cooler than LO is associated with activity (including aurorae) driven by 
the breakdown of corotation like they occur on a smaller scale for Jupiter, while the stellar radio 
emission may be a signature of high-energy particles accelerated in that process. A dimensional 
analysis, developed by Fan et al. (1982), applied to the case of L-type stars supports this hypothesis 
quantitatively 

I thank G.L. Siscoe, J.L. Linsky, and B. de Pontieu for helpful discussions and for their com- 
ments on the manuscript. 
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Table 1: Characteristic stellar properties 



Property 


M7V 


dL2 


dL5 


Stellar radius 


R* (Rq) 


0.16 h 


0.09 c 


0.08 c 


Effective temperature 


T*,eff (K) 


2700 a 


2080 fe 


1700 6 


Mass 


M* (Mq) 


0.10 h 


0.08 c 


0.07 c 


Surface gravity 


log g 


3.9 


5.4 


5.4 


Bolometric luminosity 


log( Lh i) 

CD \ L/L/t / 


30.7 


29.7 


29.3 


Moment of inertia 


log(I) d 


51.4 


50.8 


50.6 


Equatorial rotation velocity 


v eq (km/s) 


10 6 


20 6 


30 b 


Rotational energy 


log(W rot ) 


43.0 


43.5 


43.7 


Rotation period 


P(d) 


0.8 


0.23 


0.14 


Change in rotation period 


P/P (Gyr) 


3 e 


T 


> 10 e 


Synchronous rotation 


Rs/R* 


10.2 


7.3 


5.7 


Power from ISM wind 


log(P«) 


18.4 


17.8 


17.7 


Relative luminosity in Ha 


\og(L Ha /L bo i) 


-4.3 b 


-5.7 6 


< -t 


Power in Ha 


\og{L Ha ) 


26.4 


24.0 


< 22.3 


Power in X-rays 


\og{L x ) 


< 27.7 a 


< 25.5 a 


< 25. l a 


Excess of radio to X-ray 


5(L R /L X ) f 


~ xlO 


~ x3000 


? 


Loss of rotation energy 


log(W rot ) 


26.0 


26.1 


< 26.2 


Est. max. loss of rotation energy 


log(Wc) 


<27.4 


<25.9 


<25.8 


Power est. from Eq. (I13j) 


log(W) 


24.9 - 25.9 


24.8 - 25.8 


25.1 - 26.1 



a From Mohanty et al. (2004); 

b Reiners and Basri (2008); West et al. (004); 

c Baraffe et al. (2003), at age of 5 Gy; 

d For an approximate gyration radius of 0.3P*, cf. Claret and Gimenez (1990); 
e Reiners and Basri (2008), their Fig. 9. 

* Deviation from radio- X-ray relationship for warmer stars; Berger et al. (2008), their Fig. 9; 
9 Audard et al. (2007); 
h from Allen (1972). 



